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1.  Introduction 


The  prediction  of  aerodynamic  coefficients  for  projectile  configurations  is  essential  in  the 
assessment  of  the  performance  of  new  designs.  Accurate  determination  of  aerodynamics  is 
critical  to  the  low-cost  development  of  new,  advanced  guided  projectiles,  rockets,  missiles,  and 
smart  munitions.  Fins,  canards,  and  jets  can  be  used  to  provide  control  for  maneuvering 
projectiles  and  missiles.  The  flow  fields  associated  with  these  control  mechanisms  for  the  Army 
weapons  are  complex,  involving  three-dimensional  (3-D)  shock-boundary  layer  interactions,  jet 
interaction  with  the  free  stream  flow,  and  highly  viscous  dominated  separated  flow  regions  {1 
through  3).  The  jet  interference  can  extend  over  significant  portions  of  the  projectile  and  must 
be  modeled  correctly.  For  missiles,  jet  thrusters  have  been  studied  over  a  number  of  years  to 
provide  high-speed  aerodynamic  control.  These  thrusters  interact  with  the  surrounding  flow 
field,  and  the  resulting  jet  interaction  flow  field  again  is  complex.  Recently,  several  studies  have 
shown  that  tiny  synthetic  unsteady  jets  can  significantly  alter  the  flow  field  and  pressure 
distributions  for  airfoils  and  cylinders  (4,5,6).  These  synthetic  jets  are  active  control  devices 
with  zero  net  mass  flux  and  are  intended  to  produce  the  desired  control  of  the  flow  field  through 
momentum  effects.  Many  parameters  such  as  jet  location,  jet  velocity,  and  actuator  frequency 
can  affect  the  flow  control  phenomenon. 

Smith  and  Glezer  (4)  have  conducted  an  excellent  study  of  the  flow  control  by  synthetic  jets  to 
provide  increased  fundamental  understanding  of  the  flow  physics.  Amitay  et  al.  (5) 
experimentally  investigated  flow  separation  control  on  a  cylinder  using  synthetic  jet  actuators. 
Their  work  showed  that  the  interaction  of  the  synthetic  jet  with  the  free  stream  flow  resulted  in  a 
virtual  modification  of  the  body  shape  and  significantly  increased  the  lift  force  as  a  result  of  the 
flow  reattachment.  Aerodynamic  flow  control  over  an  unconventional  airfoil  has  also  been 
demonstrated  by  Amitay  et  al.  (6)  to  enhance  post-stall  performance  via  actuators  operating  at 
frequencies  higher  than  the  characteristic  frequency  of  the  airfoil.  The  synthetic  jets  are  also 
being  investigated  for  possible  applications  to  improve  heat  transfer,  reduce  drag,  and  enhance 
mixing  (7)  in  combustors,  etc.  The  present  analysis  involves  these  synthetic  jets  for  projectile 
aerodynamic  control.  The  emphasis  in  the  present  research  is  to  provide  insight  into  the 
interaction  of  these  unsteady  jets  with  the  free  stream  flow  and  to  determine  the  feasibility  of 
these  jets  for  aerodynamic  control  of  a  subsonic  spinning  projectile. 

Computational  and  experimental  data  for  these  jet  interactions  are  very  limited.  Simple  theories 
cannot  predict  the  complex  flow  fields  associated  with  the  jet  interaction,  and  experimental  tests 
are  very  expensive.  To  help  reduce  experimental  costs,  computational  fluid  dynamics  (CFD)  are 
being  used  to  predict  these  complex  flows  and  provide  detailed  pressure,  force,  and  moment 
data.  There  have  been  several  recent  numerical  studies  (8,9)  of  flow  separation  control  using 
synthetic  jet  actuators.  He  and  Krai  (8)  have  used  Reynolds-averaged  Navier-Stokes  (RANS) 


1 


models  to  study  the  effect  of  jet  location  and  jet-forcing  frequency  on  the  lift  and  drag  forces  on 
an  airfoil.  The  jet  actuator  was  shown  to  increase  the  time-averaged  lift  and  also  to  increase  the 
amplitude  of  oscillation.  Lee  and  Goldstein  (P)  have  used  the  direct  numerical  simulation  (DNS) 
to  model  a  two-dimensional  (2-D)  synthetic  jet.  Although  this  simulation  produced  very  good 
results,  the  application  of  DNS  to  practical  3-D  flows  of  interest  is  prohibitive  because  of  its 
computing  resources  requirement.  Even  large-eddy  simulations  (LES)  {10)  where  large  eddies 
are  computed  directly  and  the  small  scales  are  modeled  require  large  computational  cost 
compared  to  RANS  simulations.  While  the  RANS  method  works  well  for  steady  flows,  the 
accuracy  of  this  method  for  unsteady  flows  may  be  less  than  desired.  Since  the  large  energy- 
containing  eddies  are  computed  by  the  EES  method,  this  technique  is  more  capable  of  handling 
unsteady  shear  layers  and  wakes,  etc.  Recently,  therefore,  a  hybrid  approach  {11,12)  that 
combines  RANS  and  EES  has  been  developed  to  solve  practical  problems  of  interest  involving 
unsteady  flows  at  reasonable  computational  cost.  Both  RANS  and  hybrid  RANS-EES  models 
have  been  used  in  the  present  study. 

The  advanced  CED  capability  used  here  solves  the  Navier-Stokes  equations  {13  through  16)  and 
incorporates  unsteady  boundary  conditions  for  simulation  of  the  synthetic  jets  {17,18).  Also,  a 
hybrid  RANS-EES  turbulence  model  {12)  was  used  for  accurate  numerical  prediction  of 
unsteady  jet  flows.  Sahu  {19)  used  these  advanced  techniques  and  performed  numerical  flow 
field  computations  for  steady  and  unsteady  jets  for  a  non-spinning  projectile  configuration  at  a 
low  subsonic  speed.  Computed  lift  forces  attributable  to  the  unsteady  synthetic  jets  were  found 
to  match  the  experimental  data  well  {17,18).  The  present  numerical  study  is  a  direct  extension  of 
that  research  to  a  spinning  projectile  at  subsonic  speeds  and  atmospheric  flight  conditions.  The 
following  sections  describe  the  numerical  procedure  and  the  computed  results  obtained.  Results 
obtained  for  the  spinning  projectile  configuration  are  presented  at  Mach  0.24,  angle  of  attack 
0  degree,  and  a  spin  rate  of  67  Hz. 


2.  Solution  Technique 


2.1  CFD++  Flow  Solver 

The  complete  set  of  3-D  time-dependent  Navier-Stokes  equations  is  solved  in  a  time-accurate 
manner  for  simulations  of  unsteady  jets.  A  commercially  available  code,  CFD++  {14,15,16),  is 
used  for  the  time-accurate  unsteady  CED  simulations.  The  basic  numerical  framework  in  the 
code  contains  unified  grid,  unified  physics,  and  unified  computing  features.  The  user  is  referred 
to  these  references  for  details  of  the  basic  numerical  framework.  Here,  only  a  brief  synopsis  of 
this  framework  and  methodology  is  given. 

The  3-D,  time-dependent  RANS  equations  are  solved  by  the  following  finite  volume  method: 
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(1) 


—\^dV  +  ^\¥-G\dA  =  \mV 

dt  y  y 

in  which  W  is  the  vector  of  conservative  variables,  F  and  G  are  the  inviscid  and  viscous  flux 
vectors,  respectively,  H  is  the  vector  of  source  terms,  V  is  the  cell  volume,  and  A  is  the  surface 
area  of  the  cell  face. 

For  low  speed  flows  considered  here,  “preconditioned  implicit  relaxation”  scheme  is  used  to 
achieve  faster  convergence.  It  combines  three  basic  ideas:  (1)  implicit  local  time  stepping,  (2) 
relaxation,  and  (3)  preconditioning.  Preconditioning  the  equations  ideally  equalizes  the 
eigenvalues  of  the  inviscid  flux  Jacobians  and  removes  the  stiffness  arising  from  large 
discrepancies  between  the  flow  and  sound  velocities  at  low  speeds.  The  use  of  an  implicit 
scheme  circumvents  the  stringent  stability  limits  suffered  by  their  explicit  counterparts,  and 
successive  relaxation  allows  cells  to  be  revised  as  information  becomes  available  and  thus  aids 
convergence.  These  features  of  the  code  have  been  extremely  useful  in  the  present  numerical 
simulations  at  very  low  subsonic  speeds.  Second  order  discretization  was  used  for  the  flow 
variables  and  the  turbulent  viscosity  equation.  The  turbulence  closure  is  based  on  topology- 
parameter-free  formulations.  Two-equation  and  higher  order  hybrid  RANS-LES  turbulence 
models  were  used  for  the  computation  of  turbulent  flows.  These  models  are  ideally  suited  to 
unstructured  bookkeeping  and  massively  parallel  processing  because  of  their  independence  from 
constraints  related  to  the  placement  of  boundaries  and/or  zonal  interfaces. 

For  computations  of  unsteady  jet  interaction  flow  fields  that  are  of  interest  here,  dual  time 
stepping  (as  described  next)  was  used  to  achieve  the  desired  time  accuracy.  In  addition,  special 
jet  boundary  conditions  were  developed  and  used  for  numerical  modeling  of  synthetic  jets.  Grid 
was  actually  moved  to  take  into  account  the  spinning  motion  of  the  projectile. 

2,2  Dual  Time  Stepping 

The  “dual  time-stepping  mode”  of  the  code  was  used  to  perform  the  transient  flow  simulations. 
The  term  “dual  time  step”  implies  the  use  of  two  time  steps.  The  first  is  an  “outer”  or  global 
(and  physical)  time  step  that  corresponds  to  the  time  discretization  of  the  physical  time  variation 
term.  This  time  step  can  be  chosen  directly  by  the  user  and  is  typically  set  to  a  value  to  represent 
1/100  of  the  period  of  oscillation  expected  or  forced  in  the  transient  flow.  This  time  step  is 
applied  to  every  cell  (not  separately  varying). 

An  artificial  or  “inner”  or  “local”  time  variation  term  is  added  to  the  basic  physical  equations. 
This  time  step  and  corresponding  “inner  iteration”  strategy  is  chosen  to  help  satisfy  the  physical 
transient  equations  to  the  desired  degree.  If  the  inner  iterations  converge,  then  the  outer  physical 
transient  equations  (their  discretization)  are  satisfied  exactly;  otherwise,  approximately.  For  the 
inner  iterations,  the  time  step  is  allowed  to  vary  spatially.  Also,  relaxation  with  multi-grid 
(algebraic)  acceleration  is  employed  to  reduce  the  residues  of  the  physical  transient  equations.  It 
is  found  that  an  order  of  magnitude  reduction  in  the  residues  is  usually  sufficient  to  produce  a 
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good  transient  iteration.  This  may  require  a  few  internal  iterations  to  aehieve  (between  3  and  10, 
depending  on  the  magnitude  of  the  outer  time  step),  the  nature  of  the  problem,  the  nature  of  the 
boundary  eonditions,  and  the  eonsisteney  of  the  mesh  with  respeet  to  the  physies  at  hand. 

2.3  Unsteady  Jet  Boundary  Condition 

A  large  number  of  boundary  eonditions  (BC)  are  available  and  ean  be  speeified  at  the 
appropriate  boundaries.  Eaeh  boundary  eondition  is  eneoded  as  a  basie  form  along  with  a 
eolleetion  of  modifiers.  One  partieular  BC  used  for  the  simulations  presented  here  is  a 
“oseillating  jet”  BC.  In  its  basie  form,  it  is  a  steady  inflow/outflow  BC  wherein  the  user  supplies 
the  veloeity  normal  to  the  boundary  along  with  statie  temperature  and  any  turbulenee  quantities. 
When  the  veloeity  provided  is  negative,  it  is  eonsidered  to  be  an  inflow  and  when  it  is  positive,  it 
is  treated  as  an  outflow.  In  the  ease  of  inflow,  the  statie  temperature  and  turbulenee  quantities 
are  used  along  with  the  inflow  veloeity.  In  the  ease  of  outflow,  only  the  veloeity  is  used.  At 
inflow,  the  tangential  eomponent  of  veloeity  is  set  to  zero,  and  at  outflow,  the  tangential 
eomponent  is  extrapolated  from  the  interior.  At  outflow,  all  primitive  variables  exeept  normal 
veloeity  are  extrapolated  from  the  interior.  At  inflow,  the  statie  pressure  is  taken  from  the 
interior. 

The  first  modifier  available  for  this  BC  allows  the  veloeity  to  oseillate.  The  base  veloeity  is 
multiplied  by  an  amplitude  whieh  varies  as  sin  (27cft)  where  f  is  the  frequeney  of  the  oseillation. 
Thus,  the  oseillating  veloeity  can  cycle  from  being  positive  to  being  negative  and  back  within 
each  period  (or  from  being  negative  to  positive  and  back,  based  on  the  sign  of  the  input  for  the 
basic  BC  formulation). 

A  second  modifier  permits  the  steady  or  oscillating  inflow/outflow  to  be  “on”  over  certain  time 
intervals  and  “off’  during  other  intervals.  During  “on”  periods,  the  basic  or  the  basic  multiplied 
by  the  oscillating  amplitude  multiplier  (first  modifier)  is  used.  The  user  provides  the  ranges  of 
time  during  which  the  jet  is  on.  The  user  also  provides  a  repetition  time  period  (e.g.,  the  time 
period  corresponding  to  one  spin  rotation  of  the  projectile).  Within  each  time  period,  therefore, 
there  are  sets  of  starting  and  ending  times  that  define  when  the  jet  is  on.  During  “off’  periods, 
the  amplitude  is  set  to  zero.  In  parts  of  the  cycle  when  the  jet  is  off,  the  boundary  condition  thus 
reverts  to  the  condition  of  inviscid  surface  tangency.  This  allows  the  flow  to  slip  past  the 
boundary,  as  would  exist  (in  the  form  of  a  shear  layer)  if  the  jet  were  emanating  from  a 
cavity/hole. 

2.4  Grid  Movement 

CFD++  code  has  many  features  related  to  grid  movement.  Grid  velocity  can  be  assigned  to  each 
mesh  point.  This  general  capability  can  be  tailored  for  many  specific  situations.  For  example, 
the  grid  point  velocities  can  be  specified  to  correspond  to  a  spinning  projectile.  In  this  case,  the 
grid  speeds  are  assigned  as  if  the  grid  were  attached  to  the  projectile  and  spinning  with  it. 
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A  proper  treatment  of  grid  motion  requires  eareful  attention  to  the  details  of  the  implementation 
of  the  algorithm  applied  to  every  mesh  point  and  mesh  cell  so  that  no  spurious  numerical  effects 
are  created.  For  example,  a  required  consistency  condition  is  that  free  stream  uniform  flow  be 
preserved  for  arbitrary  meshes  and  arbitrary  mesh  velocities.  This  condition  is  satisfied  in  the 
code. 

Another  important  aspect  concerns  the  way  that  boundary  conditions  are  affected  by  grid 
velocities.  Two  significant  classes  of  boundary  conditions  (slip  or  no  slip  at  a  wall,  and  far  field 
boundary)  are  considered  here.  Both  are  treated  in  a  manner  that  works  seamlessly  with  or 
without  mesh  velocities.  In  both  cases,  consider  the  contra- variant  velocity,  which  includes  the 
effect  of  grid  motion.  Effectively,  it  is  the  normal  component  of  the  velocity  relative  to  the 
mesh.  At  the  body  surface,  it  is  set  to  zero.  For  a  no-slip  wall,  one  adds  the  requirement  that  the 
tangential  component  of  the  velocity  is  equal  to  the  mesh  velocity  by  default.  Effectively,  this 
assures  that  the  velocity  of  the  mesh  is  equal  to  the  velocity  of  the  flow  at  the  body.  At  a  far- 
field  boundary,  the  sign  of  the  contra-variant  velocity  determines  inflow  (negative  sign)  or 
outflow  (positive  sign).  The  magnitude  of  the  contra- variant  velocity  is  compared  with  the  local 
speed  of  sound,  helping  to  define  one  of  four  possibilities:  supersonic  inflow,  supersonic 
outflow,  subsonic  inflow,  subsonic  outflow.  The  characteristics  theory  is  used  to  determine  what 
and  how  much  information  is  applied  as  the  boundary  condition  for  each  type.  A  careful 
treatment  of  the  grid  motion  results  in  its  use  in  a  transparent  manner. 

2,5  Hybrid  RANS-LES  Model 

Currently,  the  two  most  popular  forms  of  turbulence  closure,  namely,  ensemble-averaged  models 
(typically  based  on  the  RANS  equations)  and  EES  with  a  sub-grid-scale  model,  both  face  a 
number  of  unresolved  difficulties.  Specifically,  existing  EES  models  have  met  with  problems 
related  to  the  accurate  resolution  of  the  near-wall  turbulent  stresses.  In  the  near-wall  region,  the 
foundations  of  large-eddy  simulation  are  less  secure,  since  the  sizes  of  the  (anisotropic)  near-wall 
eddies  approach  those  of  the  Kolmogorov  scale,  requiring  a  mesh  resolution  approaching  that  of 
a  direct  numerical  simulation.  On  the  other  hand,  existing  ensemble-averaged  turbulence  models 
are  limited  by  their  empirical  calibration.  Their  representation  of  small-scale  flow  physics 
cannot  be  improved  by  mesh  refinement,  and  over  short  time  scales,  they  tend  to  be  overly 
dissipative  with  respect  to  perturbations  around  the  mean,  often  suppressing  unsteady  motion 
altogether. 

While  EES  is  an  increasingly  powerful  tool  for  unsteady  turbulent  flow  prediction,  it  is  still 
prohibitively  expensive.  To  bring  EES  closer  to  becoming  a  design  tool,  a  hybrid  RANS-LES 
approach  based  on  limited  numerical  scales  (LNS)  has  been  recently  developed  by  Metacomp 
Technologies  {12).  This  approach  combines  the  best  features  of  RANS  and  EES  in  a  single 
modeling  framework.  The  LNS  model  is  formulated  from  an  algebraic  or  differential  Reynolds- 
stress  model,  in  which  the  sub-grid  stresses  are  limited  by  the  numerically  computed  local  length 
scale  and  velocity  scale  products.  LNS  thus  behaves  like  its  parent  RANS  model  on  RANS-type 
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grids  but  reverts  to  an  anisotropic  LES  sub-grid  model  as  the  mesh  is  refined  loeally,  thereby 
reaching  the  correct  (DNS)  fine  grid  limit.  Locally  embedded  regions  of  LES  may  be  achieved 
automatically  through  local  grid  refinement,  while  the  superior  near-wall  stress  predictions  of  the 
RANS  model  are  preserved,  removing  the  need  of  ad  hoc,  topography-parameter-based  wall 
damping. 

The  LNS  hybrid  formulation  is  well  suited  to  the  simulation  of  unsteady  fiows,  including  mixing 
flows,  and  contains  no  additional  empirical  constants  beyond  those  appearing  in  the  original 
RANS  and  LES  sub-grid  models.  With  this  method,  a  regular  RANS-type  grid  is  used  except  in 
isolated  flow  regions  where  denser,  LES-type  mesh  is  used  to  resolve  critical  unsteady  flow 
features.  The  hybrid  model  transitions  smoothly  between  an  LES  calculation  and  a  cubic  k-8 
model,  depending  on  grid  fineness.  A  somewhat  finer  grid  was  placed  around  the  body  and  near 
the  jet,  but  the  rest  of  the  flow  field  was  occupied  by  a  coarser,  RANS-like  mesh. 

To  date,  the  LNS  technique  has  been  used  successfully  on  a  number  of  unsteady  flows. 

Examples  include  flows  over  cavities,  flows  around  blunt  bodies,  flows  around  airfoils  and 
wings  at  high  angle  of  attack,  separation  suppression  using  synthetic  jets,  forced  and  natural 
convection  flows  in  a  room,  and  mixing  flows  in  nozzles. 


3.  Results  and  Discussion 


Time-accurate  unsteady  numerical  computations  using  viscous  Navier-Stokes  methods  were 
performed  to  predict  the  flow  field  and  aerodynamic  coefficients  on  the  subsonic  projectile  for 
jet-on  conditions.  Three-dimensional  numerical  computations  have  been  performed  for  the 
projectile  configuration  with  jet  interaction  with  CED-l-l-  code  at  a  subsonic  Mach  number  of 
0.24,  angle  of  attack  0  degree,  and  a  spin  rate  of  67  Hz.  The  preconditioned  version  of  the 
CED-l-l-  code  was  used  to  obtain  an  efficient  numerical  solution  at  low  speed,  for  modeling  of 
the  unsteady  synthetic  jets,  both  unsteady  RANS  and  a  hybrid  RANS-LES  approach  (72)  were 
used,  for  computations  of  these  unsteady  jets,  full  3-D  computations  are  performed  and  no 
symmetry  was  used. 

The  subsonic  projectile  is  a  1.8-caliber  ogive-cylinder  configuration  (see  figure  1).  Here,  the 
primary  interest  is  in  the  development  and  application  of  Cf  D  techniques  for  accurate  simulation 
of  projectile  flow  field  in  the  presence  of  unsteady  jets.  The  first  step  here  was  to  obtain  the 
steady  state  results  for  the  same  projectile  without  the  jet.  The  converged  jet-off  steady  state 
solution  was  then  used  as  the  starting  condition  for  the  computation  of  time-accurate  unsteady 
flow  field  for  the  projectile  with  synthetic  jets.  Computations  were  also  performed  for  the  steady 
jet  cases.  The  jet  locations  on  the  projectile  are  shown  in  figure  2.  The  jet  conditions  were 
specified  at  the  exit  of  the  jet  for  both  steady  (fixed  jet  velocity)  and  unsteady  (sinusoidal 
variation  in  jet  velocity)  jets.  The  jet  conditions  specified  include  the  jet  pressure,  density,  and 
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velocity  components.  The  flow  field  inside  the  tiny  jet  cavity  is  not  computed.  For  the  unsteady 
jet  case,  time-dependent  jet  boundary  conditions  are  applied  at  the  jet  exit.  Numerical 
computations  have  been  made  for  these  jet  cases  at  a  Mach  number,  M  =  0.24,  an  angle  of  attack, 
a  =  0  degree,  and  a  spin  rate  of  67  Hz.  The  jet  width  was  0.32  mm,  the  jet  slot  half-angle  was 
18  degrees,  and  the  peak  jet  velocities  used  were  31  and  69  m/s  operating  at  a  frequency  of 
1000  Hz. 


Figure  1.  Projectile  geometry. 


Figure  2.  Aft-end  geometry  showing  the  jet  location. 

A  computational  grid  expanded  near  the  vicinity  of  the  projectile  is  shown  in  figure  3.  Grid 
points  are  clustered  near  the  jet  as  well  as  the  boundary  layer  regions  to  capture  the  high  gradient 
flow  regions.  The  computational  grid  has  211  points  in  the  streamwise  direction,  241  in  the 
circumferential  direction,  and  80  in  the  normal  direction.  The  unsteady  simulation  took 
thousands  of  hours  of  central  processing  unit  time  on  a  Silicon  Graphics  (SGI)  Origin  computer 
running  with  16  to  32  processors. 
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Figure  3.  Computational  grid  near  the  projectile. 


Unsteady  time-accurate  CFD  computations  require  huge  computer  resources.  All  the  results 
presented  here  were  obtained  from  unsteady  numerical  computations  with  a  single  synthetic  jet 
on  a  40-mm  subsonic  grenade  (figure  1).  Looking  from  the  front  (see  figure  4),  the  projectile 
spins  clockwise  at  a  rate  of  67  Hz.  The  jet  actuation  corresponds  to  one-fourth  of  the  spin  cycle 
from  -45  to  +45  degrees  with  0  degree  being  the  positive  y-axis.  The  jet  is  off  during  the 
remaining  three-fourths  of  the  spin  cycle.  Time-accurate  unsteady  simulation  even  with  a  single 
synthetic  jet  requires  a  large  amount  of  resources.  The  unsteady  CFD  modeling  technique 
required  about  600  time  steps  to  resolve  a  full  spin  cycle.  The  unsteady  synthetic  jet  operates  at 
a  high  frequency  of  1000  Hz.  For  the  part  of  the  spin  cycle  when  the  jet  is  on,  the  jet  operated 
for  approximately  for  four  cycles.  Time-accurate  CFD  modeling  of  each  jet  cycle  required  more 
than  40  time  steps.  The  actual  computing  time  for  one  full  spin  cycle  of  the  projectile  was  about 
50  hours  using  16  processors  (i.e.,  800  processor-hours)  on  an  SGI  Origin  3000  or  an  IBM  SP3 
system  for  a  mesh  size  about  4,000,000  grid  points.  Multiple  spin  cycles  and  thus,  a  large 
number  of  synthetic  jet  operations  were  required  to  reach  the  desired  periodic  time-accurate 
unsteady  result.  As  will  be  seen  later,  some  cases  were  run  for  as  many  as  60  spin  cycles, 
requiring  more  than  48,000  processor  hours  of  computer  time.  Many  flow  field  solutions  were 
saved  at  regular  intermittent  time  intervals  to  produce  “movies”  to  gain  insight  into  the  physical 
phenomenon  resulting  from  the  synthetic  jet  interactions.  Computed  particle  traces  emanating 
from  the  jet  into  the  wake  are  shown  in  figure  5  at  a  given  instant  in  time  for  M  =  0.24  and 
a  =  0  degree.  These  traces  are  colored  by  the  velocity  magnitude.  The  particle  traces  emanating 
from  the  jet  interact  with  the  wake  flow,  making  it  highly  unsteady.  It  shows  the  flow  in  the  base 
region  to  be  asymmetric  because  of  the  interaction  of  the  unsteady  jet.  The  shear  layer  from  the 
free  stream  flow  resulting  from  the  step  corner  up  stream  from  the  base  interacts  with  the 
unsteady  jet  and  deteriorates  just  a  short  distance  down  stream  from  the  jet.  The  unsteady  jet 
substantially  alters  the  flow  field  near  the  jet  and  the  base  region,  which  in  turn  affected  the 
forces  and  moments  even  at  0  degree  angle  of  attack. 
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Figure  4.  Schematic  showing  the  jet  actuation  in  one  spin  cycle  (view 
from  the  front  or  nose). 


Figure  5.  Computed  particle  traces  colored  by  velocity,  jet  on,  M  =  0.24, 
a  =  0  degree. 

The  computed  surface  pressures  from  the  unsteady  flow  fields  are  integrated  to  obtain  the 
aerodynamic  forces  and  moments.  Computed  results  have  been  obtained  with  unsteady  RANS 
(URANS)  as  well  as  the  hybrid  RANS-LES  approach  referred  here  as  the  LNS  for  the  jet-on 
conditions.  The  unsteady  jet  is  applied  while  the  projectile  is  spinning.  For  comparison,  the  jet- 
off  computations  have  been  performed  in  a  time-accurate  mode  with  the  LNS  model.  Also,  the 
jet-off  URANS  calculations  were  first  obtained  and  the  jets  were  activated,  beginning  at  time, 
t  =  28  ms.  The  computed  normal  or  lift  force  (Fy)  and  side  force  (Fz)  are  shown  in  figures  6  and 
7  as  a  function  of  time  for  two  different  jet  velocities,  31  and  69  m/s,  respectively.  These 
computed  results  clearly  indicate  the  unsteady  nature  of  the  flow  field.  These  forces  are  found  to 
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change  as  a  function  of  time.  Changes  in  the  side  foree  are  slightly  smaller  than  changes  in  the 
lift  force.  As  shown  in  figure  6,  URANS  calculations  clearly  show  the  effect  of  the  jet  on  the  lift 
and  side  forces.  When  the  jet  is  turned  off,  the  levels  of  these  forces  drop  to  the  same  levels 
before  the  jet  aetivation  corresponding  to  the  jet-off  wake  flow.  Similar  behavior  can  be  seen  in 
figure  7  for  the  bigger  jet  ease.  The  effect  of  the  jet  is  stronger,  as  evidenced  by  the  larger  peak 
magnitudes  seen  in  the  forces. 


Figure  6.  Computed  lift  and  side  forces,  jet  on,  URANS,  M  =  0.24,  jet  velocity 
=  31m/s,  a  =  0  degree. 


Figure  7.  Computed  lift  and  side  forces,  jet  on,  URANS,  M  =  0.24,  jet  velocity 
=  69  m/s,  a  =  0  degree. 


10 


Computed  results  were  also  obtained  with  the  hybrid  RANS-LES  or  LNS  approach  for  the  jet-on 
conditions.  Figure  8  shows  the  lift  and  side  forces  obtained  with  LNS  as  a  function  of  time  for 
the  bigger  jet  case.  Here  again,  the  jet  was  turned  on  at  time,  t  =  28  ms,  starting  with  an  URANS 
solution.  The  mean  values  of  both  Fy  and  Fz  forces  are  non-zero  because  of  the  jet  interaction. 
The  LNS  model  predicts  large  amplitudes  of  oscillations  in  the  forces.  For  comparison  purposes, 
jet-off  calculations  were  obtained  with  the  LNS  model,  starting  with  the  same  initial  URANS 
solution.  A  comparison  of  the  computed  lift  force  with  LNS  for  jet-off  and  jet-on  conditions  (see 
figure  9)  clearly  indicates  that  the  wake  is  unsteady.  Also,  the  amplitudes  of  oscillations  are 
large  for  jet-off  and  jet-on  cases.  Figure  10  shows  the  comparison  of  the  predicted  lift  force  with 
URANS  and  LNS  models.  As  indicated  earlier,  the  URANS  approach  clearly  shows  when  the 
jet  is  on  and  when  it  is  off  during  the  spin  cycle.  The  effect  caused  by  the  jet  for  the  LNS  case  is 
not  as  easily  seen.  It  is  hidden  in  these  oscillations.  However,  the  mean  value  of  the  lift  force 
seems  to  be  close  to  zero  when  the  jet  is  off  during  the  spin  cycle.  In  general,  the  levels  of  the 
lift  force  oscillations  predicted  by  the  LNS  models  are  larger  than  those  predicted  by  the 
URANS. 


Figure  8.  Computed  lift  and  side  forces,  jet-on,  LNS,  M  =  0.24,  jet  velocity  = 
69  m/s,  a  =  0  degree. 
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Figure  9.  Comparison  of  computed  lift  force,  jet-off  and  jet-on  (jet  velocity  = 
69  m/s),  LNS,  M  =  0.24,  a  =  0  degree. 


Figure  10.  Comparison  of  computed  lift  force,  URANS  and  LNS,  M  =  0.24,  jet 
velocity  =  69  m/s,  a  =  0  degree. 
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Based  on  the  eomparisons  for  the  non-spinning  oases,  the  LNS  method  predioted  the  level  of 
foroes  that  olosely  matohed  the  data  {19).  In  the  ourrent  study,  the  addition  of  spin  as  well  as  the 
jet  aotuation  for  part  of  the  spin  oyole  further  oomplioates  the  analysis  of  the  CFD  results  with 
LNS.  Whether  the  jet  is  on  or  off,  the  level  of  osoillations  seen  is  quite  large  and  the  effeot  of  the 
jet  oannot  be  easily  seen  in  the  instantaneous  time  histories  of  the  unsteady  foroes  and  moments. 
To  get  the  net  effeot  of  the  jet,  unsteady  oomputations  were  run  for  many  spin  oyoles  of  the 
projeotile  both  with  and  without  the  synthetio  jets.  For  the  synthetio  jet  oases,  two  jet  velooities 
of  31  and  69  m/s  (peak  values)  were  used.  The  CFD  results  are  plotted  over  only  one  spin  oyole; 
eaoh  subsequent  spin  oyole  was  superimposed  and  a  time-averaged  result  is  then  obtained  over 
one  spin  oyole.  In  all  these  oases,  the  jet  is  on  for  one-fourth  of  the  spin  oyole  (time,  t  =  0  to 
3.73  ms)  and  is  off  for  the  remainder  (three-fourths)  of  the  spin  oyole.  Figures  1 1  through  15 
show  the  time-averaged  results  over  a  spin  oyole.  Figure  1 1  shows  the  averaged  oomputed  lift 
foroe  for  the  jet-off  ease  after  5,  10,  15,  and  20  spin  oyoles,  and  as  expeoted,  it  is  olose  to  zero, 
espeoially  after  20  spin  oyoles. 

Figure  12  shows  the  oomputed  lift  foroe  again  averaged  over  many  spin  oyoles  for  the  peak  jet 
velooity  of  31  m/s.  The  jet  effeot  oan  olearly  be  seen  when  the  jet  is  on  (t  =  0  to  3.73  ms)  even 
after  5  or  10  spin  oyoles.  The  net  lift  is  -0.07  newton  beoause  of  the  jet  aotuation  and  seems  to 
have  oonverged  after  20  spin  oyoles.  For  the  remainder  of  the  spin  oyole,  the  jet  is  off;  however, 
the  effeot  of  the  jet  on  the  wake  still  persists,  and  this  figure  shows  that  lift  foroe  is  still  available 
beoause  of  the  jet  interaotion.  The  magnitude  of  the  lift  foroe  is,  of  oourse,  reduoed  when  the  jet 
is  off.  The  foroe  and  moment  results  for  a  bigger  jet  (peak  jet  velooity  of  69  m/s)  are  shown  in 
figures  13  and  14,  respeotively.  As  shown  in  figure  13,  the  lift  foroe  requires  at  least  five  spin 
oyoles  before  reaohing  a  more  oonsistent  level.  The  lift  foroe  is  on  the  order  of  0. 17  when  the  jet 
is  on  (t  =  0  to  about  4  ms),  and  then  deoreases  gradually.  The  mean  level  of  the  lift  foroe  when 
the  jet  is  off  during  the  spin  oyole  is  on  the  order  of  0.07  newton.  The  variation  in  the  pitohing 
moment  (see  figure  14)  is  very  similar  and  seems  to  settle  after  20  spin  oyoles.  Figure  15  shows 
the  oomputed  averaged  lift  foroe  after  50  and  60  spin  oyoles  for  jet  velooities  31  and  69  m/s, 
respeotively.  It  olearly  shows  that  the  larger  jet  produoes  larger  lift  foroe  than  the  smaller  jet 
when  the  jet  is  aotivated.  Again,  one  oan  see  a  oontribution  of  the  jet  to  the  lift  foroe  for  both 
oases,  even  when  the  jet  is  off  (time  from  4  to  15  ms).  The  lift  foroe  oan  be  integrated  over  time 
to  obtain  the  impulse.  Figure  16  shows  the  impulse  obtained  from  the  lift  foroe  as  a  funotion  of 
the  spin  oyoles  for  both  jets.  As  seen  here,  in  both  oases,  it  takes  about  30  to  40  spin  oyoles 
before  the  impulse  aohieves  an  asymptote  at  a  fixed  value.  The  impulse  is  larger  for  the  larger 
jet.  Figure  17  shows  the  impulse  obtained  from  both  the  lift  foroe  and  the  side  foroe  for  the 
larger  jet  ease.  The  impulse  from  the  side  foroe  drops  olose  to  zero  at  about  30  spin  oyoles,  then 
rises  somewhat,  and  reaohes  an  asymptote  at  a  value  slightly  higher  than  zero. 
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Figure  11.  Computed  lift  force  over  many  spin  cycles,  jet-off,  LNS,  M  =  0.24, 
a  =  0  degree,  spin  =  67  Flz. 
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Figure  12.  Computed  lift  force  over  many  spin  cycles,  LNS,  jet  velocity  =  31  m/s, 
M  =  0.24,  a  =  0  degree,  spin  =  67  Flz. 
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Figure  13.  Computed  lift  foree  over  many  spin  eyeles,  LNS,  jet  veloeity  =  69  m/s, 
M  =  0.24,  a  =  0  degree,  spin  =  67  Flz. 


Figure  14.  Computed  pitehing  moment  over  many  spin  eyeles,  LNS,  jet 
veloeity  =69  m/s,  M  =  0.24,  a  =  0  degree,  spin  =  67  Flz. 
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Figure  15.  Computed  lift  force  over  many  spin  cycles  for  different  jet  velocities, 
LNS,  M  =  0.24,  a  =  0  degree,  spin  =  67  Flz. 
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Figure  16.  Computed  impulse  from  the  lift  force  as  a  function  of  spin  cycles  for 
different  jet  velocities,  LNS,  M  =  0.24,  a  =  0  degree,  spin  =  67  Flz. 
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Figure  17.  Computed  impulse  as  a  funetion  of  spin  cycles,  jet  velocity  = 

69  m/s,  LNS,  M  =  0.24,  a  =  0  degree,  spin  =  67  Flz. 

The  present  computed  results  are  currently  being  used  in  six-degree-freedom  (6-DOF)  analyses 
to  determine  the  feasibility  of  the  synthetic  jets  to  provide  control  authority  for  the  projectile;  the 
early  indications  are  that  they  can.  Also,  a  firing  test  was  just  recently  completed  to  show  that 
the  synthetic  jets  can  provide  the  control  authority  needed  for  maneuvering  a  subsonic  projectile 
to  its  target. 


4.  Summary  and  Conclusions 


This  report  describes  a  computational  study  undertaken  to  determine  the  aerodynamic  effect  of 
tiny  synthetic  jets  as  a  means  to  provide  the  control  authority  needed  to  maneuver  a  spinning 
projectile  at  low  subsonic  speeds.  Computed  results  have  been  obtained  at  a  low  subsonic  speed, 
M  =  0.24  and  0  degree  angle  of  attack,  for  a  subsonic  spinning  projectile  via  time-accurate 
Navier-Stokes  computational  technique  and  advanced  turbulence  models.  The  unsteady  jet  is 
shown  in  this  and  earlier  studies  to  substantially  alter  the  flow  field  near  the  jet  and  the  base 
region  of  the  projectile,  which  in  turn  affects  the  forces  and  moments  even  at  0  degree  angle  of 
attack.  It  is  easier  to  see  the  effect  of  the  jet  in  the  RANS  predictions  for  both  jets  (jet  velocities 
of  31  and  69  m/s).  Also,  RANS  predictions  clearly  show  that  the  jet  effect  disappears  for  the 
part  of  the  spin  cycle  when  the  jet  is  off.  In  other  words,  it  reverts  to  the  usual  wake  without  the 
jet.  The  level  of  the  lift  force  predicted  by  the  RANS-LES  model  is  larger  than  that  predicted  by 
the  URANS  technique.  This  was  also  observed  earlier  in  the  unsteady  calculations  for  the  non¬ 
spinning  projectile  case  {19). 
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In  general,  it  is  more  diffieult  to  easily  see  the  effeet  of  the  jet  in  the  unsteady  RANS-LES 
results.  Time  averaging  was  required  with  these  unsteady  results  over  many  spin  eyeles.  The  net 
time-averaged  results  obtained  over  the  time  period  eorresponding  to  one  spin  eyele  elearly  show 
the  effeet  of  the  synthetie  jets  on  the  lift  as  well  as  the  side  forees.  It  was  found  that  it  takes  10  to 
20  full  spin  eyeles  before  the  results  eonverge  to  a  periodie  state.  The  effeet  of  the  jet  in  this 
ease  is  elearly  seen  when  the  jet  is  on  during  the  spin  eyele.  However,  these  results  show  that 
there  is  an  effeet  on  the  lift  foree  (although  redueed)  for  the  remainder  of  the  spin  eyele  when  the 
jet  is  off.  This  is  a  result  of  the  wake  effeets  that  persist  from  one  spin  eyele  to  another.  The 
large  jet  is  found  to  provide  larger  lift  foree  when  the  jet  is  on.  The  level  of  lift  foree  for  the  jet- 
off  period  is  almost  similar.  The  impulse  obtained  from  the  predieted  forees  for  both  jets  seems 
to  reaeh  an  asymptote  after  30  spin  eyeles  (the  impulse  for  the  larger  jet  being  larger). 

This  work  demonstrates  a  method  to  aeeurately  prediet  the  time-aeeurate  unsteady  aerodynamies 
of  synthetie  jets  on  an  M203  grenade-launehed  projeetile.  This  eapability  has  provided 
fundamental  understanding  of  fluid  dynamie  meehanisms  assoeiated  with  the  interaetion  of  the 
unsteady  synthetie  jets  and  the  projeetile  flow  fields.  The  predieted  forees  and  moments  are  now 
being  used  in  the  6-DOF  simulations  to  assess  if  these  jets  ean  provide  adequate  eontrol  authority 
to  steer  a  projeetile  to  its  target. 
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